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ABSTRACT 

In this work we investigate the cooling channels of diffuse gas (i.e. < 0.1 cm~ 3 ) in 
cosmology. We aim to identify the wavelengths where most of the energy is radiated 
in the form of emission lines or continuum radiation, and the main elements and ions 
responsible for the emission. We use a subset of cosmological, hydrodynamical runs 
from the OWLS project to calculate the emission of diffuse gas and its evolution with 
time. We find that at z = (z = 2) about 70 (80) per cent of the energy emitted 
by diffuse gas is carried by emission lines, with the continuum radiation contributing 
the remainder. Hydrogen lines in the Lyman series are the primary contributors to 
the line emission, with a share of 16 (20) per cent. Oxygen lines are the main metal 
contributors at high redshift, while silicon, carbon and iron lines are strongest at low 
redshift, when the contributions of AGB stars and supernova la explosions to the metal 
budget become important and when there is more hot gas. The ionic species carrying 
the most energy are OlII, ClI, Cm, Sill, Si III, Fell and S III. The great majority of 
energy is emitted in the UV band (A = 100 — 4000 A), both as continuum radiation and 
line emission. With almost no exception, all the strongest lines fall in this band. At high 
energies, continuum radiation is dominant (e.g., 80 per cent in the X-ray band), while 
lines contribute progressively more at lower energies. While the results do depend on 
the details of the numerical implementation of the physical processes modeled in the 
simulations, the comparison of results from different simulations demonstrates that 
the variations are overall small, and that the conclusions are fairly robust. Given the 
overwhelming importance of UV emission for the cooling of diffuse gas, it is desirable 
to build instruments dedicated to the detection and characterisation of diffuse UV 
emission. 
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1 INTRODUCTION 

The evolution of the Universe is determined by the com- 
plex interplay of a large number of processes that in- 
volve the exchange of mass and energy. A global view of 
these mass and energy flows has been tak en to investigate 
the universal content of baryonic mass JPersic fc Saluccil 
1 19921 : IFukugita et"afl Il998l ; IFukugita fc Peebles! |2004), en- 
ergy jFukueita fc Peebles! 120041 1. me tals dPettinil |2004 



IBoucM et al. 20071 ') and star formation |Madau et al.lll996f ). 
These quantities give us useful insights in to the evolution 
of the physical Universe as a whole and they help us under- 
stand how individual processes contribute to the evolution. 

In this work, we aim to investigate a key part of this 
global energy interchange: radiative cooling of diffuse gas 
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on cosmological scales. In particular, we aim to identify the 
cooling channels that dominate the energy budget of the dif- 
fuse gas that traces the large-scale structure of the Universe 
and that feeds galaxies and fuels star formation through ac- 
cretion processes. For our purposes, we define this diffuse 
matter as the sum of three main components: i) the inter- 
galactic medium (IGM), defined as the low-density gas that 
is not gravitationally bound to galaxies and their haloes; ii) 
the circum-galactic medium (CGM), which is denser than 
the IGM, resides in the immediate surroundings of galax- 
ies, and is bound to galaxy haloes; and iii) the intra-cluster 
medium (ICM), which is the hot, shock- heated gas that ac- 
cumulates in the gravitational potentials of groups and clus- 
ters. Both the CGM and the ICM are gravitationally bound 
to haloes and represent the densest fraction of diffuse gas. 
While the distinctions between these components are often 
not clear in observations, they can be more cleanly defined 
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in simulations. In the rest of the paper, we will consider gas 
in simulations as "diffuse" if its hydrogen number density is 
nn < 0.1 cm -3 . 

Thermal and kinetic energy are injected in to the diffuse 
gas by gravity, by galactic winds powered by stellar winds 
and supernova explosions, and by feedback due to black hole 
accretion (AGN feedback), cosmic rays and other processes. 
The gas cools via a variety of mechanisms and emits radia- 
tion in a broad range of wavelengths. The cooling radiation 
can in turn heat up gas locally or far away, can be absorbed 
and re-emitted at lower frequencies by dust, or (if the mean 
free path of photons is very large) can contribute to the 
background radiation field. 

Cooling processes within galaxies are thought to con- 
tribute substantially to the diffuse infrared, optical, UV 
and X-ray backgrounds and are an essential ingredient 
in galaxy formation. These are still not well understood, 
but can be somewhat constrained and quantified by ob- 
servation. The radiation produced by the cooling of dif- 
fuse gas, on the other hand, is very difficult to observe 



directly because of it s very low s urface brightness (e.g . 
Furlanetto et all 12001 12004. 120051; lYoshikawa et al l 1200 Si; 



Fang e t al. 200 51; iBertone et al.ll2010allbll201lUGoerdt et al 



201C ; iFaucher-Giguere et al l l20ld : | van de Voort fe Schavd 
20121 ; Ivan de Voordl2012f l. Currently, the IGM can there- 



fore be more effectively studied in absorption rather than in 
emission. 

In this work, we aim to provide a global view of how 
cooling occurs in diffuse gas. To this end, we use a subset of 
runs from the Ove rwhelmingly Larg e Simulations (OWLS 
hereafter) project (|Schave et al.ll201fj| ) to address a number 
of questions concerning the emission and cooling processes 
of diffuse gas, including: i) Which transitions, ions, and gas 
phases are responsible for most of the cooling radiation, and 
should thus be the focus of theoretical modelling efforts? ii) 
Which wavebands and emission lines carry most of the en- 
ergy emitted by diffuse gas through cooling and fluorescence 
radiation, and should thus be primary targets for observa- 
tional studies? iii) What is the relative importance of con- 
tinuum emission and line emission? iv) What is the global 
emission of diffuse gas, and how does it change with time? 
How much energy is stored in the IGM? Since our focus is 
diffuse gas, we will ignore emission by dust and by gas with 
densities exceeding tih =0.1 cm -3 . 

The paper is organized as follows. Section [2] describes 
the emissivity tables used to calculate the gas emission and 
shows the dependence of the emissivity on temperature for 
a large sample of ions, individual emission lines and the con- 
tinuum. Section discusses the simulations and the method- 
ology used to calculate the gas emission. Section [4] presents 
results for the fraction of energy emitted by different ions, 
elements and the continuum in the reference simulation set 
(REF in the following) at z — and z = 2. In Section [5] we 
compare results from different simulations to understand the 
effect of varying the physical prescriptions. Finally, we draw 
our conclusions in Section[6]and present a series of numerical 
convergence tests in the Appendix. 



Table 1. Adopted solar abundanc es, f rom I Al lende P rieto et all 
fcOOlI' l. lAfiende Prieto et alj <2002l) and iHolwegerl <200lft . 

Element Ui/nn Element ni/njj 



H 


1 


Mg 


3.47X1CT 5 


He 


0.1 


Si 


3.47X1CT 5 


C 


2.46X1CT 4 


s 


1.86xl0" 5 


N 


8.51X1CT 5 


Ca 


2.29xl0" 6 


O 


4.90X1CT 4 


Fe 


2.82xl0" 5 


Ne 


l.ooxicr 4 







2 THE EMISSIVITY TABLES 

In this Section we describe the emissivity tables used to cal- 
culate the contribution of different species to the total emis- 
sion. In particular, in Section \2. II we show the temperature 
dependence of the emissivity for all ions, and for a sample of 
the strongest individual emission lines for each ion; Section 
12.21 describes the emissivity tables for the continuum. 
The tables have been created using 

CLOUD-vO 

as a 

function of gas temperature (in bins of ALog 10 T = 0.05 for 
10 2 K < T < 10 8,5 K), hydrogen number density (in bins of 
ALog 10 n H = 0.2 for 10~ 8 cm" 3 < n H < 10 cm" 3 , but note 



that we will ignore emission from gas with tih > 10 cm ), 
and redshift. 

The adopted solar abundance is Zq = 0.0127, corre- 
sponding to the value obtained using the default abundance 
set of CL OUDY (see Table [I). This s e t combines the abun- 
dance s of I Allende Prieto et all l|200ll '). lAllende Prieto et all 
(|2002h and lHolweged (l200ll) and can differ significantly from 
the estimates of ILoddersI (I2003T) . In part icular, the abun- 
dance of oxygen given by ILoddersI (|2003h is about 20 per 
cent higher than in CLOUDY. Note that variations in the 
adopted solar abundances may introduce variations in the 
computed emissivity per unit gas mass, but not for the sim- 
ulation results, which are instead calculated using the abso- 
lute abundances predicted by the simulations. 

The tables created with CLOUDY assume the gas to 
be dust-free, optically thin and in (photo-)ionisation equi- 
librium in the presence of a uniform, evolving cosmic UV/X- 
ray back ground radiation contri buted by both galaxies and 
quasars (jHaardt fe Madaull200ll ) and the cosmic microwave 
background spectrum. While dust emission is an extremely 
important cooling channel for the interstellar medium, its 
contribution to the cooling of the diffuse IGM is probably 
much smaller, particularly at higher temperatures where the 
dust is efficiently destroyed. On average, local sources of ion- 
ising radiation are expected to become more important than 
the background for gas c louds with ne utral hydrogen column 
densities > 10 cm" " 2 i|Schavdl2006h , which are expected 
to have volum e densities slightly below nu — 10" 1 cm" 3 
l|Schavell200~lh , the highest gas density we consider. As ion- 
ising radiation sup presses cooling rates due to the remova l 
of bound electrons ([Efstathioull 19921 ; IWiersma et~"ai"1l2009ar ) . 
we may have slightly overestimated the cooling rates for 
high-density gas, part icularly for temperatures < 10 5 K (see 
W iersma et alj|2009ar ) . Ionisation equilibrium is expected to 



1 Version c07.02.02, last described in lFerland et al.lll998l and re- 
leased in July 2008. See http://www.nublado.org/. 
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hold for photo-ionis ed regions and for dense gas in the cen- 
tres of clusters (see IBertone. Schave fc Dolaa 20081 for a re- 
view). The gas is prevalently photo- ionised at low density 
and low temperatures and is collisionally ionised at high 
densities and temperatures. The equilibrium assumption can 
break down for shock-heated gas. This might happen, for ex- 
ample, for warm- hot gas in the outskir t s of groups and clus- 
ters, a s argued by Yoshida et al. (I2005T) : Yoshikawa fc Sasakil 
fcOOd k lCen fc Fand (|2006l ): lGnat fc Sternberg! (|2007l . l2009h . 
However, since these studies ignored photo-ionisation, they 
may have overestimated the importance of non-equilibrium 
ionisation. 



2.1 Line emissivity 

We use the line emissivity tabl es created by 



Bertone et al.l 

d2010ef>. which were also used inrBertone et all (|2010bl . l201lh 
and Ivan de Voort fc Schave! (|2012h . They include a rela- 
tively small selection of the strongest lines for each ion - 
slightly over 2000 emission lines, out of a total of several mil- 
lions - for the 11 elements explicitly tracked by OWLS (H, 
He, C, N, O, Ne, Mg, Si, S, Ca and Fe). The contributions 
of all other elements are small, and should not significantly 
affect our estimate of the total emission. For comparison, 
calcium is the element in our sample that contributes the 
smallest fraction of the total energy, of the order of 0.1 per 
cent (see Section l4~Tj) . The contributions of t he neglected ele- 
ments are all smaller than this value (see also lWiersma et alj 
l2009aD . 

We restrict our attention to the lines with the highest 
emissivities for each ion, so that the total emissivity is ac- 
curate to within about a per cent or better. In many cases 
this requires just a few lines per ion (C iv 1548A, 155lA and 
O vi 1032A, 1038A are good examples, carrying almost 100 
per cent of the total energy emitted by the C iv and O vi 
ions respectively), while for heavier atoms such as iron and 
calcium many more lines are required. 

The emissivity of the strongest lines of each ion is shown 
as a function of temperature in Figs. Q] (H, He, C, N, O, 
Ne, and Mg) and [2] (Si, S, Ca, and Fe). All emissivities are 
shown at z = for a constant hydrogen number density 
of «h = 10" 3 cm" 3 and solar element abundances. Note 
that emissivity tables are thus computed for a metallicity 
that substantially exceeds the average meta llicity of diffuse 
gas, w hich is about 0.1 Zq at low redshift (|Wiersma et al.l 
l2009bl ). For densities lower than the «h = 10 _3 cm -3 as- 
sumed here, some of the emissivity curves show higher val- 
ues in the low temperature regi me, where photo-ionisa tion 
becomes dominant, as shown by IBertone et all l |2010al ). At 
higher densities photo-ionisation becomes negligible, and the 
emissivity curves drop sharply for temperatures below the 
values corresponding to the peak emissivities, as is also the 
case for most metal lines for the density shown here. For 
clarity, for each ionic species we show at most the three 
lines with the highest peak emissivities, even if more lines 
are included in the tables. When less than three lines are 
shown, as for the UV doublets of lithium-like ions and for 
beryllium-like ions, the lines that are not shown have peak 
emissivities that are at least an order of magnitude smaller 
than the dominant lines. In some cases, as for most iron and 
calcium lines, we show only one line for each ion to avoid 
overcrowding the figures. 



The total line emissivities of individual ions, calculated 
by summing up all the lines included in the tables for each 
ion, are shown in Fig. [3]as a function of temperature, at red- 
shift z = 0, for «h = 10~ 3 cm" 3 and solar abundances. The 
colour of the emissivity curves of each ion varies smoothly 
with increasing ionisation state; black-blue colours show the 
emissivities of lower ionisation states with low peak tem- 
peratures, while yellow-red colours represent highly ionised 
species, such as helium- and hydrogen-like ions. For these 
last two species, the emissivity peak roughly indicates the 
temperature threshold above which the atoms are nearly 
fully ionised. With the exception of iron lines, the colour- 
coding used for Figs. [T][2] and Fig. [3] is the same. 

Fig. [3] illustrates visually how the shape and the peak 
temperature of the line emissivity curves depend on the 
atomic number and on the ionisation state of each element. 
In particular, the peak temperature, i.e. the temperature for 
which the emissivity is maximum, increases steadily with in- 
creasing ionisation state and atomic number. This can easily 
be understood, because more bound electrons require larger 
energies, and therefore higher temperatures, to be removed 
or excited. Hydrogen and helium contain at most two elec- 
trons on the Is energy level that can be extracted at rela- 
tively low temperatures. Heavy elements with atomic num- 
bers up to 11 (C, N, O, and Ne) have two electrons on the 2s 
level and up to six electrons on the 2p level; this translates 
into emissivity curve peaks that cluster in two separate tem- 
perature ranges: those for helium- and hydrogen-like ions at 
the highest temperatures, and all others for ions with re- 
maining electrons on the 2s or 2p levels clustering at lower 
temperatures. Elements with atomic numbers between 11 
and 18 (which includes Mg, Si, and S) contain electrons on 
the 3s and 3p levels. These produce a further set of emis- 
sivity curves that cluster together and peak at the lowest 
temperatures. Finally, calcium and iron have electrons on 
the 4s and 3d levels respectively and for this reason a fourth 
and a fifth separate cluster of emissivity curves, correspond- 
ing to the upper energy levels, appear at the lowest temper- 
atures. Moving up the scale in atomic mass, each new set 
of emissivity curves, corresponding to a higher energy level, 
always appears at the lowest temperature, while the whole 
distribution moves to higher temperatures, with lines from 
hydrogen-like ions having the highest peak temperatures. 

As we will sh ow in the following S ections, and as al- 
ready discussed by IBertone et all (|2010al lbT). the shapes and 
peak temperatures of the emissivity curves are fundamen- 
tal in determining the level of the line emission. However, 
the relative intensity of the emissivity curves does not auto- 
matically translate into equal relative emission intensities, 
because the total emission also depends on the ionic abun- 
dance of different species and on the overall mass content 
in each temperature range. As an example, the emissivity 
of Lya lines is lower than those of many oxygen lines, but 
overall the total emission from hydrogen lines is much larger 
than that of oxygen lines, because hydrogen is much more 
abundant than oxygen. This will be shown in Section 0] 

2.2 Continuum emissivity 

In this Section we describe our technique to build the con- 
tinuum emissivity tables, which is similar to that used for 
the line emissivity tables. 
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z=0.00 n H =10" 3 cm" 



-26 



3 -28 



-22 



-24 




H 1 1216A 
H 1 1026A 
H 1 972.5A 



.4- 




. C 1 1656A 

C 1 1 561 A 

. C 2 1335A 

C 2 892.8A 

C 2 2328A 

. C 3 977.0A 

C 3 1907A 



-4- 




C 4 1548A 

C 4 1551A 

C 5 41.47A 

C 5 40.27A 

C 5 40.73A 

C 6 33.74A 

C 6 28.4 7A 



I 

N 5 1239A 

N 5 1243A 

N 6 29.53A 

N 6 28.79A 

N 6 29.08A 

N 7 24.78A 

N 7 20.9 1A 



3 "28 



-22 




He 1 625.6A 
He 1 584.3A 
He 1 537.0A 
He 2 303.8A 
He 2 256.3A 
He 2 1640A 



26 ' ~ 




. N 1 1200A 

. N 2 907.1A 

N 2 1085A 

N 2 634.4A 

. N 3 991.0A 

N 3 678.2A 

, N 4 765.0A 

N 4 I486A 




-24 



-26 



5 -28 



-22 




. 1 6300A 

I 6363A 

. 2 833.8A 

2 425.2A 

2 530.8A 

. 3 698.2A 

3 500.7A 

3 835.0A 




4 549.3A 

4 789.0A 

5 630.0A 

5 1218A 

6 1032A 

6 1038A 



-24 



-26 



3 -28 



-22 




. Ne 2 461. IA 
Ne 2 444.2A 
Ne 2 12.81/1 

. Ne 3 486.9A 
Ne 3 3869A 
Ne 3 3968A 

. Ne 4 546.2A 
Ne 4 2424A 




Ne 5 573.4A 
Ne 5 481.2A 
Ne 6 400.0A 
Ne 6 561.8A 

Ne 7 464.9A 
Ne 7 97.77A 




Ne 8 770.4A 
Ne 8 780.3A 
Ne 9 13.45A 
Ne 9 13.70A 
Ne 9 13.55A 
NelO 12.14A 
NelO 10.24A 



-24 



-26 



3 -28 




. Mg 1 2853A 
. Mg 2 2796A 
Mg 2 2803A 

. Mg 3 183.5A 
. Mg 4 322.7A ' 
. Mg 5 3S3.2A 
. Mg 6 406.4A 



. Mg 7 276.9A 
. Mg 8 315.5A 

. Mg 9 368.5A 
. Mg10 609.8A 
Mg10 624.9A' 



. Mg11 9.169A 
Mgll 9.314A 
Mgll 9.231A 

. Mg12 8.422A 
Mgl2 7.106A 




5 6 7 8 
Log, T (K) 



5 6 7 8 
Log, T (K) 



5 6 7 8 
Log, T (K) 



Figure 1. The emissivity of the strongest individual emission lines for each ion as a function of temperature for njj = 10 cm - , 
z = 0, and solar abundances. Lines from the same ion are displayed with the same colour. Lines are shown for hydrogen, helium, carbon, 
nitrogen, oxygen, neon and magnesium. 



The total continuum emissivity can be separated into 
the contribution from gas purely composed of hydrogen 
and helium, and contributions from each of the metals 
tracked in OWLS. The emissivity tables for heavy elements 
are computed u sing C LOUDY, following the approach of 
W iersma et al.l (|2009ar ). The emissivity due to metal y is 



taken to be the difference between the emissivity computed 
assuming solar abundances and the emissivity computed af- 
ter setting the abundance of y to zero, while keeping the 
abundances of all oth er elements the s ame re lative to hydro- 
gen. As discussed bv lWiersma et all (|2009al ) . this approach 
is accurate as long as the contribution of heavy elements 
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z=0.00 n H =10" 3 cm - 



"I I I I" 



-22 



-24 



-26 



3 -28 




. S! 2 34.81m 
Si 2 1814A 

. St 3 1207A 
S! 3 1883A 

. Si 4 1 394A 
Si 4 1403A 

. Si 5 95.55A 
. Si 6 249.1A 



.4- 





-22 



-24 



-26 



3 -28 




2 750.3A 

3 667.0A 

4 655. 1A 

4 737.3A 

5 786.5A 

5 1198A 

6 9J3.4A 

7 59.83A 

8 203.1 A 



-22 



-24 



-26 



3 -28 




. Co 1 4228A 
. Co 3 351.3A 
. Co 4 3I0.4A 

. Co 5 314.5A 
. Co 6 330.4A 
. Co 7 338.0A 
. Co 8 356.7A 



Co 9 470.7A 
Co 10 576.0A 
Coll 30.42A 

Co 12 28.08A 
Co 1 3 169.4A 
Co14 24.15A 




. Co 15 22.85A ' 
. Co16 171.7A" 

. Co 17 204.2A . 
. Co 18 302.0A ■ 
177A" 
. Co20 3.023A ' 




-22 



-24 



-26 



3 -28 




. Fe 2 2400A 
Fe 2 25.99m 
. Fe 3 1126A 

. Fe 4 532.9A 
. Fe 5 3892A 
. Fe 6 279.2A 
. Fe 7 163.5A 



-22 



-24 



-26 



3 -28 




Fe13 204.3A 
Fe14 219.8A 
Fe15 284.2A 

Fe16 335.4A 
Fe17 15.26A 
Fe17 17.05A 
Fe17 16.88A 



5 6 7 8 
Log, T (K) 



5 6 7 8 
Log, T (K) 



5 6 7 8 
Log 10 T (K) 



Figure 2. As Fig. [T] but for silicon, sulphur, calcium and iron. 



to the free electron density is small, i.e. for metallicities 
Z < Zq. Since hydrogen and helium are both important 
contributors of free electrons, we tabulate the total emis- 
sivity for these elements combined and do so for a range of 
helium abundances (relative to hydrogen). 

Hence, while the line emissivity tables are 3- 
dimensional, the continuum tables are effectively 4- 



dimensional for metals and 5-dimensional for the combined 
hydrogen and helium continuum. As for the lines, the contin- 
uum emissivity is sampled as a function of density, temper- 
ature and redshift. One additional dimension is due to the 
wavelength dependence of the continuum, which is mapped 
over 12 orders of magnitude in wavelength in the range 
A = 0.01 A to 1 m. The tables for the contribution of hy- 
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z = 0.00 



n H =10" 5 cm" 3 




4 5 6 7 
Log, T (K) 



4 5 6 7 
Log, T (K) 



4 5 6 7 8 
Log, T (K) 



4 5 6 7 
Log, T (K) 



Figure 3. The total line emissivity of ions as a function of temperature for njj = 10 cm - , 2 = 0, and solar metallicity. Each panel 
shows ions of a specific element. The line colours vary as violet-blue-green-yellow-orange-red with increasing ionic number. The dotted 
red line shows the total line emissivity of each element. The bottom-right panel shows the total line emissivities of each element. The 
temperature for which the emissivity of each line is maximum increases with ionic number. With the exception of iron, the colour-coding 
of ions is the same as in Figs. [T]- [2] 



drogen and helium include one more dimension, the 5th, 
that takes into account variations in the abundance of he- 
lium relative to that of hydrogen. Tables have been created 
for 8 different helium mass fractions, varying in the range 
Y = 0.238 - 0.308 by bins of AY = 0.01, with Y = 0.248 
being the primordial abundance. 

The key factor determining the level and the shape of 
the continuum emissivity is the gas temperature. In Fig. [3] 
we show how the total continuum emissivity for a gas with 
solar composition varies as a function of temperature be- 
tween T = 10 3 K and T = 10 s K. The overall shape of the 
curves reflects the wavelength and temperature dependence 
of the emissivity for thermal Bremssstrahlung (i.e. free-free 
emission), which scales as T~ 1//2 exp(— hc/XkT), where h is 
Planck's constant, c is the speed of light and k is Boltz- 
mann's constant. For T < 10 6 K the emissivity curve be- 
comes less smooth at intermediate wavelengths as radiative 
recombinations and two-photon decays of metastable levels 
become increasingly important for lower temperatures (e.g. 



Osterbrocklll974l ; ISutherland fc Dopitdll993l ; iFerland et all 
Il998h . From left-to-right (i.e. from short to long wave- 
lengths), the sharp peaks in the 10 3 K curve correspond to 
the Hen Lyman limit (228A, n = 1), the Hei Lyman limit 
(504A, n = 1), and the Hi Lyman (912A, n = 1), Balmer 
(3646A, n = 2), and Paschen (8204A, n = 3) limits. The 
sharp peaks at longer wavelengths correspond to recombi- 
nations to higher energy levels (n > 3). The broader peak 
at « 2431 A corresponds to 2-photon emission from Hi. 



We do not show the contributions of the individual ele- 
ments to the continuum emissivity as we find that it is nearly 
always dominated by hydrogen and helium. The exception 
is the regime hc/X S> kT; however in this regime the con- 
tinuum emissivity is exponentially suppressed. Finally, we 
note that dust emission, which we have ignored here as it is 
not clear how much dust resides in the IGM, is expected to 
become increasingly important at longer wavelengths. 
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mass weighted metallicity 




-8 -6 -4 -2 -8 -6 -4 -2 

Log n H (cm -3 ) Log n H (cm"') 



Figure 5. The mass-weighted distribution of gas in the temperature-density plane (left panel, contours are spaced by 1 dex) and the 
mass-weighted average gas metallicity (right panel) for the REF_L100N512 run at z = 0. The vertical, dashed lines indicate the star 
formation threshold, above which we impose an effective equation of state onto the gas. The gas metallicity is not a good tracer of the 
gas mass, but is highly biased towards high densities. 
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Figure 4. The total continuum emissivity as a function of tem- 
perature for constant density, njr = 10 — 3 cm -3 , 2 = 0, and solar 
abundances. The overall shapes of the curves reflect the wave- 
length and temperature dependence of thermal Bremsstrahlung. 
The peaks at intermediate wavelengths that become prominent 
for low temperatures correspond to radiative recombinations and 
2-photon emission from hydrogen and helium. 



3 SIMULATIONS 

In this Section we briefly describe the suite of cosmologi- 
cal, hyd rodynamical simula tions that constitute the OWLS 
project (jSchave et al.ll20ld ). with particular attention to the 
model referred to as the "reference", or REF, model. 

The OWLS project aims to investigate the effects of 



varying a large number of physical prescriptions on the pre- 
dicted evolution of the baryons in the Universe. At present, 
more than fifty simulations have been run, with varying mass 
resolutions, box sizes, and physical prescriptions for gas cool- 
ing, star formation, chemical evolution, supernova and AGN 
feedback. The simulations have been run with a custom ver- 
sion of GADGET III, a significantly extended version of the 
parallel PMTree Smoothed Par ticle Hydrodynamics (SPH) 
code GADGET II |Springej|2005l ). 

Our basic results have been computed using two simu- 
lations of the reference model, which contain 2 x 512 3 par- 
ticles each. The first, labelled REF '_L025iV512, has a box 
size of 25 /i" 1 comoving Mpc and follows the evolution of 
structures down to z = 2. The gas and dark matter parti- 
cle masses ar(0 1.4 x 10 6 h' 1 M and 6.3 x 10 6 h' 1 Mq, 
respectively, and the gravitational force resolution (i.e. the 
softening) is 0.5 h~ x comoving kpc, but is frozen at 0.5 h~ x 
proper kpc below z = 2.91. The second simulation, denoted 
by REFJ.1QQN512, has a box size of 100 h' 1 comoving Mpc 
and follows the evolution of structures down to z = 0. Com- 
pared with REF_L025N512, this simulation uses 4 3 = 64 
times larger particle masses and a 4 times larger softening. 
We have tested additional runs with variations in the physi- 
cal prescriptions, box sizes and resolution, which we describe 
in more detail in Section [5] and in the Appendix. 

Besides gravitational and hydrodynamical forces, the 
simulations employ new physical prescriptions for radia- 
tive cooling, star formation, chemodynamics, supernova and 
AGN feedback. Radiat ive cooling is impleme nted follow- 
ing the prescription of IWiersma et all (|2009al ). The cool- 
ing rates are calculated for each element under the assump- 
tion of an optically thin gas in ionisation equilibrium, in 



2 These are the initial values. The masses of baryonic particles 
vary slightly during the course of the simulation as a result of 



stellar mass loss. 
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the presence of the cosmic microwave background and an 
ionising background ra diation from quasars and galaxies 
l|Haardt fc Madaull200ll ). Pre-computed tables for each ele- 
ment have been created with CLOUD\ff], as described for 
the continuum emissivity in Section \2. 2 1 The totalgas cool- 
ing rate is calculated by interpolating the tablecl and by 
summing up the contributions of each element, rescaled to 
the particle element abundance. 

Star formation is imp lemented as described in 
ISchave fc Dalla Vecchial l|200Sl ). Gas with densities nu > 
1CP 1 cm -3 is expected to be multi-phase and unstable to 
star formation. As we lack both the physics and the reso- 
lution to model the multiphase interstellar medium, we im- 
pose a polytropic equation of state onto these particles of 
the form P oc p 4 ^ 3 . This choice of power-law index results 
in a constant Jeans mass, which prevents numerical frag- 
mentation. The observed Kennicutt-Schmidt star formation 
law is analytically converted into a pressure law and gas 
particles with hh > 1CF 1 cm -3 are stochastically converted 
into star particles at th e pressure-dependen t rate implied by 
this la w. As shown by ISchave et al .1 (|2010l ) and lHaas et al.l 
(2012), the star formation rates are insensitive to the as- 
sumed star formation law as a result of self-regulation. We 
will not consider any simulations with variations in the star 
formation prescriptio ns. The prescrip t ion for chemical evo- 
lution is described in IWiersma et al l (|2009bl ). The scheme 
tracks the evolution of the abundance of the 11 elements that 
contribute significantly to the radiative cooling rate (H, He, 
C, N, O, Ne, Mg, Si, S, Ca, and Fe) produced by massive 
stars (core collapse supernovae and stellar winds) and in- 
termediate mass stars (Type la super novae and asymp totic 
giant branch [AGB] stars), assuming a lChabrierl (2003) stel- 
lar initial mass function. 

Feedback from core collapse supernovae and stellar 
winds from massive stars are included as kinetic feedbac k 
following the prescription of lDalla Vecchia fc Schavd l|200^) . 
The model describes winds using two free parameters, the 
mass- loading factor 77 = m w /m„, which parametrizes the 
rate at which mass is injected into the wind, a func- 

tion of the star formation rate m*, and the wind initial ve- 
locity, « w . The REF model assumes r\ = 2 and i> w = 600 km 
s~\ which correspond to about 40 per cent of the energy 
available from core collapse supernovae. Model REF does 
not include AGN feedback. 

The mass-weighted distribution of gas in the 
temperature-density plane is shown in the left panel 
of Fig. [5] for the REF_L100N512 run at z = 0, while the 
right panel shows the mass- weighted average gas metallicity. 
The distributions for the REF.L025N512 run (not shown) 
are qualitatively similar, but the gas overall fills a more 
compact region of the plane. In particular, the maximum 
gas temperatures at z = 2 are lower than at z = by 
about half a dex, and the minimum densities are higher by 
a similar factor. 

At both redshifts, the bulk of the metals do not trace 
the bulk of the IGM mass. Instead, most of the gas mass is 
concentrated along the temperature-density relation corre- 
sponding to the low-density, photo-ionised IGM (red contour 



3 version 05.07 rather than 07.02 as in lWiersma et al I J2009ah . 

4 Using equation (3) rather than (4) of lWiersma et at (2009a). 



in the left panel of Fig. [5]), while metal-rich gas concentrates 
in the circum-galactic medium, with nrr > 10 -5 cm -3 and 
T < 10 5 K, and in the intra-group and intra-cluster me- 
dia, with n H > 10~ 5 cm~ 3 and T > 10 6 K. A significant 
amount of metals can finally be found in the low-density, 
warm-hot IGM with 10 5 K < T < 10 6 K. A more in depth 
st udy of the metal distrib ution in these runs can be found 
in lWiersma et all |2009bl ). 

3.1 The gas emission 

T he gas line emis s ion is e stimated following the methodology 
of iBertone et al. I |2010al ). The luminosity of each particle is 
calculated by interpolating the emissivity tables as a func- 
tion of density and temperature at the desired redshift and 
is then rescaled to the corresponding element abundance. 

The luminosity Lij [z) of particle i corresponding to 
line I of element y at redshift z, in units of erg s~ , is calcu- 
lated as 

i,,,(z). £l , ( 2 ,r„«H,)^^, (i) 

where m gaSi ; is the gas mass and pi the density of particle i. 
£;,q (z,T, bh) is the line emissivity in units of erg cm -3 s -1 
and is bi-linearly interpolated (in logarithmic space) from 
the CLOUDY tables as a function of the particle tempera- 
ture T; and hydrogen number density nn,i- Xyj is the mass 
fraction of element y and X y Q its solar value. For our cal- 
culations we use the "smoothed" element abundances de- 
scribed in W lersma et alJ (|2009bl) . which are defined as the 
ratio of the SPH smoothed metal mass density and the SPH 
smoothed mass density, Z y< i = p y ,i/pi- This definition of the 
gas metallicity is more consistent with the SPH formalism 
than the commonly used particle metallicity, i.e. the ratio 
of the metal and total m ass of the particle. As discussed 
in W lersma et al l (|2009bl ). the use of smoothed abundances 
has the additional advantage that it reduces the effects of the 
lack of metal mixing that is inherent to the SPH technique. 

Similarly to the line emission, the continuum emis- 
sion is calculated by interpolating the tables to the desired 
density, temperature, and redshift. The tables for the con- 
tributions of the metals are interpolated individually and 
rescaled to the corresponding elemental abundances. The 
hydrogen+helium tables are in addition interpolated to the 
helium abundance of the particles. As a result, the contin- 
uum luminosity as a function of wavelength and redshift 
L c , y ,i (z, A) for metal y and particle i, in units of erg s 
Hz" 1 , is given by: 

L a ,y,i (z, A) = e y , A , (z,T it nn,i) , (2) 

Pi X y Q 

with ej/,A,0 (z, T, nn) the continuum emissivity of metal y in 
units of erg cm -3 s _1 Hz -1 . The continuum emission for 
hydrogen+helium is 

£c,H+Hc,i (z, A) = £H+Hc,A (z, Ti,rlH,i, Xuc,i) — - — — , (3) 

Pi 

with £h+Hc,a (z, T, riH, Ahc) the combined continuum emis- 
sivity for hydrogen and helium, also in units of erg cm -3 s _1 
Hz - . The total continuum emissivity of particle i is then 

L c ,i (z, A) = L c , H +Ho,i (z, A) + ^2 Lc,y,i (z, A) . (4) 

y 
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We will often plot the volume- weighted, specific emis- 
sion E v , which is also an average emissivity, computed by 
summing the luminosities of all gas particles in the simula- 
tion and dividing by the simulation volume. We will compute 
this specific emission for selected lines and continuum com- 
ponents, for the sum of all lines of particular ions, and for 
the sum of all emission processes from all ions. Finally, if we 
only include gas particles that fall in a certain temperature 
range, we can compute these quantities as a function of the 
gas temperature. 



4 THE COOLING CHANNELS OF DIFFUSE 
GAS 

In this Section we present results for the two reference sim- 
ulations: i?.£F JL025N512 at z = 2 and i?£F_L100N512 at 
z = 0. In Section 14.11 we discuss the relative contributions 
of ions and the continuum to the global emission budget, 
and in Section [472] we investigate the main gas cooling chan- 
nels as a function of temperature. Finally, in Section T4. 31 we 
present results for the global energy density of the diffuse 
IGM. 

4.1 Cooling channels integrated over all diffuse 
gas 

In this Section we discuss the relative contributions of dif- 
ferent ions and the continuum to the global emission budget 
of diffuse gas, integrated over the entire simulation volume. 

Fig. E] shows results for the fraction of energy emitted 
in the form of lines by each iorjf], relative to the total en- 
ergy in lines carried by the corresponding element, for all 
elements explicitly followed in the OWLS runs. The total 
line energy emitted by all ions in each panel equals one. Re- 
sults are shown at z = for _R_EF_L100N512 (red diamonds) 
and at z = 2 for ii£F_L025N512 (blue, filled circles). The 
bottom right panel of Fig. [6] shows the fraction of energy 
emitted in the form of lines by each element, with the last 
two columns showing the total fraction in lines and con- 
tinuum. This panel demonstrates that lines and continuum 
carry similar amounts of energy, with lines being slightly 
more important. Of the emision lines, those from hydrogen 
carry most of the energy, their contribution to the total emis- 
sion being about thirty per cent. At z = the next most 
important elements are silicon, carbon, helium, and iron, 
whose lines contribute ~ 3 — 6 per cent each. Next come 
oxygen, nitrogen, sulphur, magnesium, and neon which con- 
tribute w 0.6 — 2 per cent each. Finally, calcium lines only 
account for about 0.08 per cent of the energy. 

Looking now at the other panels, we see that there are 
a few common features to all elements and ions. Most of 
the energy in lines is emitted by atoms in relatively low 
ionisation states, which indicates that dense, warm gas is 
responsible for much of the energy dissipated in the form of 
line radiation. Lines from neutral gas in general contribute 

5 The emissivity tables do not include any lines for a few ionic 
species, because their cmissivities are very small. The neglected 
ions are Nei, Mgm, Si I, Siv, Si, Sxi, Sxn, Can and Fei. The 
corresponding curves are thus missing from Figs. l9l and 1101 



Table 2. Wavelength bands. 



Band 


A min (A) A max (A) 


X-ray 


0.1 


100 


Far ultraviolet (FUV) 


100.0 


1000 


Near ultraviolet (NUV) 


1000.0 


4000 


Optical 


4000.0 


7000 


Infrared and longer (IR+) 


7000.0 


1 m 



a minor fraction of the emission, with the notable exception 
of Oi, which contributes about 15 per cent of the energy 
emitted by oxygen at z = (but is unimportant at z = 2). 

At z — 2 the relative importance of lines and contin- 
uum is unchanged, but there are significant changes in the 
contributions from individual elements. Helium and oxygen 
are much more important and each account for nearly ten 
per cent of the energy. This comes at the expense of silicon, 
nitrogen, and iron which account for about 3 — 7 times less 
energy than at z = 0. The fact that helium lines contribute 
less at lower redshift can be understood by noting that the 
metallicity of the diffuse gas increases with time, providing 
a larger range of cooling channels in the relevant tempera- 
ture range. The change in relative contributions of the heavy 
elements is also partly a consequence of the chemical evolu- 
tion of the diffuse gas. While core collapse supernovae dom- 
inate the metal production rate early on, intermediate mass 
stars release large amounts of nitrogen and iron on longer 
time scales through A GB winds and type la supernovae (e.g. 
IWiersma et~alll2009rj ). 

The evolution in the relative contributions of individual 
ions is larger than that for elements. The lowest ionisation 
states are much less important at higher redshifts, prob- 
ably because there is less dense gas and because the UV 
background is more intense. The highest ionisation stages of 
oxygen and heavier elements carry much smaller fractions 
of the energy at higher redshifts. This is because at high 
redshift the fraction of hot, metal enriched gas is very small 
compared to the mass of wa rm gas (see e.g. ICen fc Ostrikerl 
1 1999 ; IWiersma et~alll2009bh . 

Fig.[7]shows how the fraction of energy emitted through 
continuum radiation and lines is distributed in the different 
wavelength bands defined in Table [2] The left panel shows 
results for AEF_L100N512 at z = and the right panel for 
i?£F_L025N512 at z = 2. 

Independent of redshift, almost half the energy is emit- 
ted in the FUV band and another third in the NUV band. 
As such, about 80 per cent of the total energy is emitted 
in the UV wavelength range A = 100 - 4000 A. The opti- 
cal emission contributes 5-8 per cent of the energy, and the 
IR+ band about 7-15 per cent. The band with the largest 
variation with redshift is the X-ray band: its contribution 
increases from less than 1 per cent at z — 2 to about 7 per 
cent at z = 0. 

The 7Ji5F_L025N512 simulation has a box size of only 
25 ft -1 Mpc and it is possible that the fraction of energy 
emitted in the X-ray band is underestimated because the 
most massive haloes are not sampled. However, since the 
fraction of hot gas with T > 10 s K does not exceed a few 
per cent in mass at z — 0, we do not expect to have under- 
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Figure 6. Fraction of the energy emitted per unit volume by (almost) all ions for each element, averaged over all the diffuse gas in the 
simulation. In each panel a different element is shown and the sum of the contributions of all ions corresponding to that element equals 
1. The bottom right panel shows the fraction of energy per unit volume emitted by each element, by all lines and by continuum, with 
respect to the total specific emission. Results are for the i?FF_L100N512 simulation at z = and for the RFF_L025N512 simulation at 

2 = 2. 



estimated the contribution of the X-ray band by more than 
a few per cent. 

Overall, the continuum emission dominates the energy 
budget at the highest energies (X-ray and FUV bands), 
while emission lines are dominant at longer wavelengths 
(A > 1000 A). The relative contributions of lines and contin- 
uum vary little with redshift. The exceptions are lines with 
A > 7000 A, whose contribution increases from 4 to 12 per 
cent of the total over time, tripling their share in the IR+ 
band. 

The evolution of the sums of the volume-averaged co- 
moving emissivities of all emission lines from individual el- 
ements, along with the total radiative cooling rates in the 
forms of lines and continuum, are shown for 2 = — 5.5 in 
Fig. H for simulation _R£F_L100N512. The global emissivity 
increases by about a factor of 3 from z = 5.5 to 2, after 
which it drops increasingly rapidly to a value that is similar 
to the z = 5.5 value. The shape of the curve is very similar to 



that of the cosmic star formation history in this model (see 
ISchave et al"1l2010r ). consistent with the idea that the rate at 
which diffuse gas cools ultimately driving the rate at which 
galaxies form stars. As for the star formation rate, the evo- 
lution is driven by structure formation, which increases the 
density contrast and thus the cooling rates (most of which 
scale as the density squared), and the general expansion of 
the Universe, which decreases the gas densities and slows 
down the formation of structure, particularly when the en- 
ergy density of the Universe ceases to be dominated by mat- 
ter. 

Results for the _R£F_L025N512 (not shown) and 
i?£F_L100N512 simulations converge well below z ~ 4. 
At higher redshifts, the 7?i?F_L025N512 simulation predicts 
significantly higher global emissivities. The high-redshift 
divergence of these two simulations, which again paral- 
lels t hat of the cosmic star formation rates ijSchave et alJ 
2010), occurs because the higher-resolution run can resolve 
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REF_L100N512 z=0 REF_L025N512 z = 2 




Figure 7. Fraction of the energy emitted in the continuum emission in different wavebands, averaged over all the diffuse gas in the 
simulation. Results are for the RFF_L025N5f2 simulation at z = 2 (left panel) and for the iJ£FX100N512 simulation at z = (right 
panel). In both cases, more than half the energy is emitted in the FUV band, and a large fraction of the remaining energy is emitted 
in the NUV band. The X-ray band becomes significant only at low redshift, where it emits more energy than the optical and 1R bands 
combined. Here, the IR band includes all wavelengths larger than 7000 A (i.e., including the proper IR and all other bands up to radio 
frequencies). 




2 4 6 8 

redshift 

Figure 8. Evolution of the specific emission (emission per unit 
volume) of lines from individual elements and from the con- 
tinuum, averaged over all the diffuse gas. Results are for the 
R&FX100N512 simulation. 



lower mass galaxies that form earlier than galaxies in the 
REF _L100N512 run and that can enrich the IGM with met- 
als at higher redshifts. The greater fraction of high-density 
gas and the higher metallicity both increase the global, high- 
redshift cooling rate. 



Continuum and line emission each account for about 
half the emission over the whole redshift range, with lines 
becoming relatively more important with time. Hydrogen 
accounts for about half of the line emission at all redshifts. 
While helium is by far the second most dominant element 
at high redshift, for z < 2 oxygen is responsible for a similar 
amount of line emission. At z < 1 the contributions of all 
elements heavier than hydrogen are within a factor of three 
of each other, except for neon, magnesium, and particularly 
calcium, which are less important in terms of the global 
emission budget. 

The relative strength of the emission from different 
elements varies with time as a consequence of chemical 
evolution, evolution in the distribution of mass in the 
temperature-density plane, and evolution in the UV back- 
ground radiation. The elements that differ the most from 
the general trend are iron and silicon, whose global emis- 
sivity continues to increase down to z ~ 0.5. For iron this 
is probably due to a combination of two facts: the increase 
of the fraction of hot gas, where iron emission is most ef- 
ficient, and the increase in iron produced by type la su- 
pernovae, which release their nucleo-synthetic products on 
much longer time scales than elements such as oxygen and 
neon that are predominantly released by massive stars. The 
production of silicon is, however, also dominated by mas- 
sive stars. To understand its evolution, we note that at low 
redshift its emission is totally dominated by the Sin 34.8/^ 
line, which is produced by cold T ~ 10 3 K gas (see Figs. \E\ 
and llip . Hence, it is the red cloud of high- metallicity, high- 
density, low-temperature gas in the right panel of Fig. [5] 
that dominates the emission. At high redshift the gas has a 
lower metallicity and less of it is able to cool down to the 
temperatures where Sin becomes important. 
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REF_L100N512 



z = 0.00 




Figure 9. Fraction of the energy emitted per unit volume in the form of lines as a function of temperature, averaged over all the diffuse 
gas in the Ri7F_L100N512 simulation at z = 0. The different curves show different ions and the different panels show different elements, 
except for the bottom-right panel, which compares the fractions of the total mean emission per unit volme due to lines (red, dashed) and 
continuum (solid, blue). The fractions are relative to the total emission in each temperature bin. The results were obtained by averaging 
the emission over the entire volume of the Ri?F_L10uN51 2 simulation at z = 0. The curves end at the maximum temperature sampled 
by the particles in the simulation volume. 



4.2 Cooling channels as a function of temperature 

In this Section we discuss the main cooling channels of dif- 
fuse gas as a function of temperature. 

Fig. [9] shows the relative contributions of the total line 
emission from individual ions to the average energy that is 
emitted per unit time and volume as a function of the gas 
temperature. The results are for 7?i5F_L100N512 at z = 
and have been averaged over all the diffuse gas in the entire 
simulation volume. Different panels correspond to different 
elements, as indicated, and different coloured curves cor- 
respond to different ions. For each ion, all lines have been 
summed. The curves are normalised separately for each tem- 
perature, i.e. the solid curves show the fraction of the total 
emission from gas at that particular temperature emitted 
by a single ion. The dashed red curve in each panel shows 
the sum of all the ions of the corresponding element, i.e., it 
is the sum of all the solid curves in the same panel. Finally, 



the dashed red and the solid blue curves in the bottom-right 
panel show the fractions of the total energy that is emitted 
in the form of lines and continuum, respectively, summed 
over all elements. The two curves therefore sum to unity. 

Fig. [10] is similar to Fig. [9] but shows the results for 
i?£F_L025N512 at z = 2. The two simulations differ in 
terms of their resolution and box size. However, as the re- 
sults are nearly converged with respect to both resolution 
and box size (see the Appendix), they can be directly com- 
pared in order to check for evolution. 

It is important to note that Figs. [9] and [10] do not sim- 
ply reflect the emissivity curves. As the emission is averaged 
over the entire simulation volume, the results depend on the 
distribution of mass and the chemical composition of the gas 
as a function of temperature and density. The figures thus 
illustrate the relative contributions of line and continuum ra- 
diation (bottom-right panel) and the relative contributions 
of different elements and ions to the cooling rates of diffuse 
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Figure 10. As Fig. [9j but for the i?EF_L025N512 simulation at z = 2. 



gas as a function of temperature after averaging over the 
entire (simulated) universe. 

Focusing first on the bottom-right panel of Fig. [9] we 
see that emission lines carry most of the energy at tem- 
peratures lower than a few times 10 6 K, while continuum 
emission is dominant for higher temperatures, where it can 
exceed the contribution of lines by more than two orders of 
magnitude. This is because most elements are very highly 
ionised at such high temperatures, thereby severely reducing 
the number of bound electrons that can be collisionally ex- 
cited. Bremsstrahlung therefore dominates the cooling rates 
of such high-temperature plasmas. Line and continuum radi- 
ation contribute approximately equally to the cooling rates 
for T ~ 10 4 K, but lines dominate for both lower and higher 
temperatures. At T ~ 10 5 K the continuum only accounts 
for about ten per cent of the radiative losses. 

Fig. shows that among emission lines, H i lines carry 
most of the energy for T < 10 4 K, while metal lines are domi- 
nant for T > 10 5 K. All elements that we consider contribute 
fractions 3> 1 per cent of the total energy in some tempera- 
ture range. The elements (ions) whose lines contribute sub- 
stantially more than ten per cent of the total energy are 



hydrogen for T < 10 4 K (H i), carbon at ~ 10 4 7 K (C n and 
Cm), oxygen for ~ 10 5 K (Om, O iv, Ov, and Ovi), neon 

10 3 K 
10 6 K 



at 10 5 b K (Nev), magnesium for 10 5,9 K, silicon at 



(Sii), w 10 45 K (Sin), and 



' 10 6 K, and iron for 
10 7 K. 



(Fe viii and Fe ix) and at T s 

Comparing Figs. [9] and 1101 we find little evolution in 
the distributions, with the main difference being the relative 
contributions of lines and the continuum in the range 10 4 ' 5 < 
T < 10 6 K. In this temperature interval, the contribution 
of lines is about twice as large at z — than at z = 2. The 
increase is due to the higher average metallicity of the IGM 
at low redshift. 

Figs. [11] and [12] show the contributions of individual 
lines, rather than ions, at z = and z — 2 respectively. All 
lines that contribute more than 1 per cent of the total emit- 
ted energy in at least one temperature bin are shown and 
they are listed in Tables [Bl] and [B2] in Appendix[B] The Ta- 
bles show the total fraction of energy radiated through each 
line along with the fraction of energy carried by the line at 
the temperature for which the line is most important. Both 
quantities have been sorted by the fraction of the energy 
that they carry. 
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Figure 11. Fraction of the energy emitted per unit volume in the form of individual lines as a function of temperature. The fractions 
are relative to the total emission in each temperature bin. Only lines that contribute at least 1 per cent of the total emitted energy in 
at least one temperature bin are shown. The results were obtained by averaging the emission of diffuse gas over the entire volume of the 
REFX100N512 simulation at z = 0. 



The lines that account for more than ten per cent of 
the total energy are Hi 1216A for T < 10 4 K, Hen 304A at 
« 10 4 9 K, Cm 977A at ~ 10 48 K, Oiv 549A at « 10 5 2 K, 
Ov 630A at » 10 5 - 4 K, O vi 1032A at 10 5 5 K, Nevi 400A at 
10 5 5 K, Sin 34.8^ at ~ 10 3 K, and Feix 170A at « 10 5 ' 9 K. 
Comparing Figs. llll and ll2l which show z — and 2, respec- 
tively, we see that the same lines dominate, but that most 
of the metal lines are less important at higher redshift and 
that the reverse is true for the Hen lines. 

Fig. Q3] shows how, for each temperature, the fraction 
of energy emitted through lines and continuum radiation is 
distributed over the wavelength bands defined in Table [2] 
The left and right panels show results for 7?EF_L100N512 
at z = and fl.EF_L025N512 at z = 2, respectively. Note 
that the sum of all wavelength bins yields the curves shown 
in the bottom-right panels of Figs. R?l and 1 1 01 

As the temperature increases, a greater fraction of the 
energy is emitted at shorter wavelengths. For T 2> 10 6 K 
the X-ray continuum is by far the most important cooling 
channel, while UV lines account for the vast majority of ra- 



diative losses for 10 4 K <C T < 10 6 K. If we sum line and 
continuum radiation, then the UV dominates down to the 
lowest temperature plotted, T = 10 3 K. IR+ lines only be- 
come important at T ~ 10 3 K, where they account for about 
30 per cent of the emitted energy. The maximum contribu- 
tion of the optical band is about 10 per cent and is reached 
near 10 4 K, with most of the energy coming out in the form 
of lines. Comparing the two panels, we see little evolution, 
although the FUV continuum is more important for gas with 
T ~ 10 5 K at high redshift. 



4.3 The global energy density of the IGM 

Finally, we estimate the importance of the cooling radia- 
tion emitted by diffuse gas to the cosmic energy budge t 
of the Universe. According to iFukugita fe Peebles! l|2004D . 
the contribution of the IGM to the cosmic energy bud- 
get due to kinetic energy associated to the gas peculiar 
velocities is of order Qigm, kinetic ~ 10 -8 . Indeed, we find 
fiiGM.kinctic ~ 5 x 10" 8 for _R_£F_L100N512 at z = and 
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Figure 12. As Fig.[TT] but for the _R_BF_L025N512 simulation at 



« 0.9 x 1(T 8 for REF _L025N512 at z = 2. These values are 
comparable to the total internal energy, for which we mea- 
sure fiiGM, thermal ~ 3 x 10 -8 at both redshifts. Low density 
gas with p/pmean < 10 carries about half the internal energy 
at z — 2 and about 10 per cent at z = 0. 

The fact that the total kinetic and internal energy den- 
sities are comparable is not surprising. In regions where the 
kinetic energy density exceeds the thermal one, peculiar ve- 
locities will typically be supersonic, resulting in the con- 
version of kinetic to thermal energy through shocks. If, on 
the other hand, the thermal energy density is substantially 
larger than the kinetic energy, e.g. as a resuls of explosive 
feedback, then the gas will tend to expand at the velocity 
of sound until this velocity becomes comparable to the local 
escape velocity, which will in turn be similar to the typical 
peculiar velocities in the region. 

By dividing the cosmic thermal energy density by the 
global cooling rate, we infer a cosmic cooling time of about 
4.1 Gyr at z — and 1.4 Gyr at z = 2. For both redshifts the 
radiative cooling times are about a factor of three smaller 
than the age of the Universe. Thus, the diffuse IGM as a 
whole would tend to radiate a large fraction of its internal 
energy over a Hubble time, if all heating sources were re- 



moved. The fact that the radiative cooling time is compara- 
ble to the expansion time scale probably reflects the fact that 
structure formation, and therefore heating through accre- 
tion shocks and feedback from star formation, proceeds on 
this time scale and that radiative cooling and shock-heating 
approximately match each other. 

Figure [T4l shows the total emitted energy per unit vol- 
ume as a function of the gas density and temperature. The 
left panel shows results for the _R_Ei ? _L100N512 simulation at 
2 = 0, the right panel for the i?_B_F_L025N512 at z = 2. The 
total emitted energy in each density and temperature bin 
has been calculated as the sum of the continuum emission, 
integrated over all wavelengths, plus the energy emitted by 
all lines. When compared to Fig. [5] Fig. [14] demonstrates 
that the bulk of the emitted energy traces neither the bulk 
of the mass nor the bulk of the metals in the simulations. In- 
stead, Fig. [Til can be bet t er und erstood by c omparing it to 
Fig. 13 of iBertone et al.1 (l2010ah. Fig 10 o f IBertone etUI 
i|2010bl ) and Fig. 9 of IBertone et al.l (1201 ll ). These figures 
show the temperature-density distributions of gas, weighted 
by the emission in different UV and X-ray lines, and indi- 
cate how the bulk of the emission is produced by the denser 
gas, independently of the gas mass distribution in the case 
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Figure 13. The mean fraction of the total energy emitted by diffuse gas per unit volume accounted for by emission lines (dashed) and 
continuum radiation (solid) and by different wavelength bands (colours; sec Table [2] for definitions). The results are averaged over the 
entire simulation volume. The left and right panels are for REFX100N512 at z = and iJi?F_L025N512 at z = 2, respectively. For 
T > 10 6 ' 5 K most of the energy is emitted as continuum radiation in the X-ray band, while UV emission lines dominate the radiative 
cooling rate at 10 4 < T < 10 6 ' 5 . For T ~ 10 3 K the sums of all emission lines falling in the IR+, FUV, and NUV bands each account for 
about 30 per cent of the radiative energy losses. There is little evolution. 



REF_L025N512 z = 2.00 REF_L100N512 z = 0.00 




Log n H (cm -3 ) Log n H (cm -3 ) 

Figure 14. The total emitted energy per unit volume as a function of the gas density and temperature. The left panel shows results for 
the i?EFX100N512 simulation at Z = 0, the right panel for the i?i?F_L025N512 at z = 2. Most of the energy is emitted by relatively 
dense and metal rich gas. 
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of the metal lime emission. Similarly, Fig.[5]shows that most 
of the total emission at both redshifts is produced by rela- 
tively dense gas, in accordance with the oc p 2 dependence in 
the cooling law. The main difference between the results at 
z — and at z — 2 is the increasing importance of emission 
from warm-hot and hot gas (T > 10 6 K) with decreasing 
redshift, as seen also in Figs. [6] and [71 



5 SIMULATIONS WITH VARYING PHYSICS 

In this Section we compare predictions from a subset of sim- 
ulations from the OWLS project with different implementa- 
tions of a number of physical proces ses. The full set of O WLS 
simulations is described in detail in ISchave et all l|20lti ) and 
we refer the interested reader to that work. 

We use two subsets of simulations: a set with box sizes 
100 ft -1 Mpc for results at z = 0, and a set with box sizes 25 
ft" 1 Mpc for results at z — 2. All runs contain 2 x 512 3 parti- 
cles and use the same particle masses and initial conditions 
as the reference runs. Wherever possible, we use simulations 
with 25 and 100 ft" 1 Mpc boxes that implement the same 
physical prescriptions at both redshifts. In the following we 
only give a brief description of the runs we consider for our 
comparison. 

(i) REF: Our reference model, as described in Section 

(ii) NOZCOOL: As the REF model, but with gas cooling 
rates calculated assuming primordial element abundances. 
The neglect of the contributions of metal lines implies slower 
gas cooling than in the REF model. Note that we do not 
turn off metal-line cooling in post processing the simulation 
to predict the emission, which is thus not self-consistent. 

(iii) NOSN: As the REF model, but without supernova 
feedback: the metals produced by stars are only transported 
by gas mixing and no energy is transferred to the IGM and 
ISM when SNe explode (100 ft" 1 Mpc box only). 

(iv) NOSN.NOZCOOL: As the REF model, but without 
supernova feedback and with gas cooling rates for primordial 
element abundances (25 ft" 1 Mpc box only). 

(v) NO A GB_NOSNIa : As the REF model, but with no 
chemical enrichment by supernova la and AGB stars and no 
supernova la energy transfer (100 ft" 1 Mpc box only). 

(vi) AGN: As the REF model, with the addition of 
A GN feedback imp l emen ted according to the prescriptions 
of iBooth fc Schavd l|2009l ). 

(vii) WVCIRC: Model with a momentum-driven imple- 
mentation of galactic winds , following the prescriptions of 
lOppenheimer fc Pavel (|2008l ). with the difference that winds 
are produced "locally" to the star formation event and are 
fully hydrodynamically coupled as in the REF model. The 
wind initial velocity and the wind mass loading factor are 
defined 3<7 and r\ — u w o/c respectively, with 
a = v2v irc the galaxy velocity dispersion, v c i IC the halo 
circular velocity, estimated using an on-the-fly friends-of- 
friends halo finder, and u w o = 150 km s" 1 . 

Fig. H5l shows the specific emission of ions, elements and 
the continuum for simulations with box sizes 100 ft" 1 Mpc 
at z = 0. Fig. 1161 shows the same quantities for simulations 
with box sizes 25 ft Mpc at z — 2. The left panels show 
the specific emission of the ions most relevant to the global 



cooling rate, while the right panel shows the contributions of 
individual elements, all lines combined and the continuum. 

Qualitatively, Figs. [15] and [16] confirm the results of 
Fig. [5] that lines carry a larger fraction of the energy ra- 
diated away by cooling processes than the continuum, and 
that most of the cooling is produced by hydrogen lines. Vari- 
ations in the relative contributions of lines and the contin- 
uum to the global budget are minimal for all models, with 
the expception of the NOSN model. Values of the specific 
emission of ions and elements vary among simulations by up 
to two orders of magnitude in the most extreme cases, with 
the largest variations for the NOZCOOL model at z = and 
for highly ionised species (e.g. Ovm). 

The NOZCOOL model implements a less efficient cool- 
ing prescription than other models. As a consequence, the 
gas temperature is artificially high, with the largest vari- 
ations at low redshift, where the gas is on average more 
metal-rich and metal-line cooling is more severely underesti- 
mated. While the relative fractions of energy in lines and the 
continuum do not show significant variations, the fraction of 
energy emitted by hydrogen and helium lines is higher in the 
NOZCOOL model than in REF, while the fractions emitted 
by metals are lower by up to a dex. When ions are con- 
sidered, the difference tends to be larger for highly ionised 
species that emit energy in the X-ray and FUV bands than 
for the lowest ionisation species. 

Results for the NOSN.NOZCOOL model are roughly 
similar to those of REF at z — 2. However, results for the 
NOSN model are significantly different from those of REF 
at z = 0. Line emission has an overall stronger impact than 
continuum emission for this model, and the individual con- 
tributions to the budget of ions and elements present signif- 
icant variations. In particular, some cooling channels (e.g. 
silicon lines) are enhanced, while others (e.g. oxygen and 
magnesium lines) are suppressed. 

The contribution of metal-line emission to the global 
budget is noticeably lower for the A GN model than for REF. 
This is because AGN feedback quenches star formation and 
suppresses the production of metals, which translates into a 
lower specific emission for metal lines. 

Variations in the emission channels due to different feed- 
back mechanisms are smaller than for other physical varia- 
tions. This is demonstrated by the results for the WVCIRC 
model, which are on average within a factor of a few of those 
of REF and present the smallest variations in our sample of 
models. This means that the way metals are injected in to 
the IGM through supernova feedback is largely irrelevant in 
the determination of the global cooling channels. 

In the NOAGB.NOSNIa model, the production of met- 
als by AGB stars and type la supernovae is turned off, which 
implies that the relative abundances of some elements, such 
as carbon and iron, are lower than in the REF model. This 
directly translates in to smaller relative contributions to the 
gas cooling by these same elements, with differences of up 
to 0.5 dex. 

In conclusion, the comparison between different simu- 
lations demonstrates that: i) The gas cooling rates are the 
most essential ingredient to determine the gas temperature 
and ultimately the relative importance of different cool- 
ing channels. Since the cooling rates depend on metallic- 
ity, which on average increases with decreasing redshift, the 
underestimation of the gas cooling rate when assuming pri- 
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Figure 15. The specific emission (emission per unit volume) at z = for a subset of OWLS runs with boxes of 100 h^ 1 Mpc on a 
side. The left panel shows the ions with the largest contributions to the total emission, while the right panel shows the contributions of 
individual elements, all lines combined and the continuum. 



z = 2 L=25 h"' Mpc 



-2 7 



o 

°, -3 



-4 7 



i 

A 

* O w 

c I 



o 



• 

5K 



• 



u 







+ 




-5 



>-'-CSICS|fO'*CS|N->'*00<NM->CM (NK1 IO(NIOCO 

<D <D (DOJCP d)0)(U 

IIIOOOZOOOZZ2tOtOLOU.lj.lj. 



O 



* 



# REF 

* NOSN_NOZCOOL 
NOZCOOL 

AGN 
WVCIRC 



§ « 

$ i $ 

O O 



+ 



o 



cn ._ o <u 

5 CO CO O Ll. 



59 »k» 



ions elements 

Figure 16. As Fig. [HI but for a subset of OWLS runs with boxes of 25 h 1 Mpc at z = 2. 



mordial abundances becomes more severe with decreasing 
redshift. ii) The largest differences between the predictions 
of different simulations are visible for highly i onised species, 
such as Oviii. This confirms the findings of iBertone et al.l 
(|2010ah that X-ray observations might provide a strong tool 
to rule out or verify the theoretical predictions from different 
models. 



6 CONCLUSIONS 

In this work we have investigated the cooling channels of dif- 
fuse gas in the Universe. We limited our analysis to cooling 
radiation from diffuse, non star forming gas with tih < OT 
cm -3 , and we neglected the emission from stars, supernova 
explosions, active galactic nuclei and compact objects. 

We used a subset of high resolution, cosmologi- 
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cal hydrodynamica l simulations from the OWLS project 
ijSchave et alj|2010h . with box sizes of 25 ft" 1 Mpc and 100 
h^ 1 Mpc on the side and 2 x 512 3 particles. Lower resolu- 
tion simulations have been used to test the dependence of 
the results on resolution and box size (see Appendix) . Most 
of our results were inferred rom the REF simulations, which 
include element-by-element cooling for 11 species, star for- 
mation, stellar evolution and chemodynamics, and galactic 
winds driven by star formation. 

Our main results can be summarised as follows: 

(i) Emission lines carry about 80 per cent (70 per cent) 
of the energy emitted by cooling radiation at redshift z — 2 
(z = 0), with the remainder carried by continuum emission. 
Of the energy emitted by lines, the majority is carried by 
Lya and other Lyman series lines of hydrogen. Other im- 
portant coolants are oxygen, helium and carbon at z = 2, 
and silicon, carbon, iron and helium at z = 0, in order of 
decreasing contributions. This shift in the relative contribu- 
tions of different elements reflects variations in the produc- 
tion of metals with time: while oxygen from massive stars is 
the most abundant metal at high redshift, the production of 
carbon and iron by AGB sta rs and supernova la ex plosions 
is boosted at lower redshifts (|Wiersma et al. I l2009bh . In ad- 
dition, structure formation increases the amount of hot gas 
with time, favouring the line emission from elements with 
higher atomic numbers, such as iron. 

(ii) The emission lines that carry most of the energy usu- 
ally correspond to relatively low ionisation states of atoms. 
In particular, some of the strongest lines correspond to O in, 
Cn, Cm, Sin, Sim, Fen and Sin ions. The wavelengths of 
the strongest lines of each ion can be read in Figs. [T] and [2] 
and the strongest lines for each element are shown in Figs. 
1111 and [T2l and listed in Appendix [Bl 

(iii) Emission in the UV band strongly dominates the en- 
ergy budget at all redshifts. With almost no exceptions, the 
strongest emission lines are UV lines. Overall, the contribu- 
tion of the continuum emission to the energy budget is larger 
than the contribution of line emission at high energies (X- 
ray and FUV bands), while lines dominate at lower energies 
(NUV, optical and IR+ bands). X-ray emission carries a 
few per cent of the total energy, its contribution increasing 
with decreasing redshift. The optical and IR+ bands to- 
gether contribute less than 20 per cent of the total emission, 
with IR+ lines dominating the budget at z = 0. 

(iv) We have compared results for a number of OWLS 
runs that implement variations of a selected number of phys- 
ical prescriptions. We find that the gas cooling rates are the 
most essential ingredient to determine the gas temperature 
and ultimately the relative importance of different cooling 
channels. Since the cooling rates depend on the gas metal- 
licity, which on average increases with decreasing redshift, 
the underestimation of the gas cooling rate when assuming 
primordial abundances becomes more severe at low redshift. 
The largest variations in the predictions preferentially affect 
highly ionised species, such as Ovm. 

Our results demonstrate that the UV band is a trea- 
sure trove of information about the cooling processes that 
determine the evolution of the diffuse IGM. Current IGM 
studies mostly target absorption lines in QSO spectra such 
as Lya, Civ, Ovi. These lines are extremely useful and are 
relatively easy to detect. However, Civ and Ovi lines only 



carry a small fraction of the total energy radiated away by 
cooling gas, compared with other UV lines from lower ioni- 
sation states, such as C in, O ill and O v, which account for 
a much larger portion of the energy budget. 

The fact that most of the cooling radiation is emitted 
in the UV band, both by the continuum and by emission 
lines, may have some interesting consequences. In fact, pho- 
tons emitted at wavelengths shorter than the Lyman limit 
break at 912 A have a very short mean free path at hig h 
redshift (e.g. Ilnoue fc Iwatal [20081 ; iProchaska et al.1 12009 ). 
As such, a large fraction of the energy originally emitted in 
the FUV band is quickly absorbed and re-radiated at longer 
wavelengths and it might be problematic to detect cooling 
radiation from diffuse gas directly. 

At low redshift, however, the mean free path of ionising 
photons increases sub stantially, and the Universe becomes 
more transparent (e.g. lMadau et al"1ll999T ). This means that 
a larger fraction of the cooling radiation in the FUV band 
could be detected. Indeed, evidence of FUV emiss ion in clus- 
ters has been detected by EUVE observations l|Lieu et al.l 
1996). It would therefore be highly desirable to build a satel- 
lite with the capability to try to detect emission at high spec- 
tral resolution in the FUV band A ~ 100 — 1000 A, where 
most of the strongest emission lines and a large fraction of 
the continuum emission are found. 
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APPENDIX A: CONVERGENCE TESTS 

In this Appendix we investigate how resolution (Appendix 
IA1[) and box size f Appendix lA2|) affect our results. We show 
results for the specific emission of oxygen ions, which are 
broadly representative of all elements in all bands. Equiva- 
lent results hold for all elements traced by OWLS. In Ap- 
pendix [S3] we show how the binning in wavelength used to 
calculate the continuum emission affects the estimate of the 
continuum contribution to the total energy budget. All re- 
sults are shown for runs that use the REF physical model. 



Al Resolution 

We investigate the effect of varying the numerical resolution 
by comparing results for the specific emission of oxygen ions 
in simulations with varying particle number and constant 
box sizes. The simulations contain N — 2 x 128 3 , 2 x 256 3 
and 2 x 512 3 particles each. As a consequence, the mass and 
the spatial resolution vary in steps of 8 and 2, respectively 
Fig. IA1I shows results for simulations with 25 h" 1 Mpc 

0' 



boxes and DM particle masses Mum = 4.1 x 10 s h~ L M 
5.1 x 10 7 h' 1 M and 6.3 x 10 6 h' 1 M at z = 2 (left 
panel) and for simulations with 100 h~ l Mpc boxes and 
DM particle masses Mom = 2.6 x 10 10 h' 1 Mq, 3.2 x 10 9 
h' 1 M and 4.1 x 10 s h' 1 M at z = (right panel). 

Results are well converged at z — 2 for all particle 
masses and deviations among different sets of simulations 
are minimal. At z — only the two higher resolution runs 
converge well, while the run with N = 2 x 128 3 particles 
diverges significantly for most ions. This is n ot surprising, 



since neither the star formation rate density (Schaye et all 
nor the line emission (e.g. iBertone et alj|2010al ) con- 
verge for this simulation. The convergence in the 25 hT 1 Mpc 
simulations at z — 2 for particle masses equal or smaller than 
those used in the 100 ft -1 Mpc run with TV = 2 x 512 3 parti- 
cles, suggests that a good level of convergence is reached in 
the highest resolution run with a box size of 100 h' 1 Mpc. 
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Figure Al. Dependence of the specific emission (emission per unit volume) on mass resolution. Results are shown for runs with L = 25 
h _1 Mpc at z = 2 (left panel) and for runs with L = 100 h —1 Mpc at z = (right panel), for constant box size. The number of 
particles in the simulations is N = 2 X 128 3 , 2 X 256 3 and 2 X 512 3 , corresponding to DM particles masses Mdm = 4.1 x 10 s h _1 Mq, 
M DM = 5.1 X 10 7 h' 1 M and Af DM = 6.3 X 10 6 h' 1 M (left panel) and to A/ DM = 2.6 X 10 10 h" 1 Mq, M dm = 3.2 X 10 9 h' 1 Mq 
and M DM = 4.1 x 10 8 h~ l Mq (right panel). 
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Figure A2. Dependence of results on the box dimension L. Re- 
sults are shown at z = for simulations with constant resolution 
and box sizes of 25, 50 and 100 ft -1 Mpc. 



A2 Box size 

In this Appendix we investigate the effect of varying the box 
dimension at constant mass resolution. We compare results 
for simulations that have DM particle masses A^dm = 4.1 x 
10 8 h' 1 Mq and box sizes of 25, 50 and 100 h' 1 Mpc at 



z — 0. The number of particles in each run is TV = 2 x 128 3 , 
N = 2 x 256 3 and TV = 2 x 512 3 , respectively. 

Results for the specific emission of oxygen are shown in 
Fig. IA2I and show good convergence for all box sizes. O iv 
emission shows variations of about 0.5 dex for the 25 ft -1 
Mpc run, also found in nitrogen and carbon ions at equiv- 
alent ionisation temperatures. Small differences are seen in 
highly ionised species of calcium and iron (e.g. Fexxv and 
Fexxvi), which emit in the hard X-ray band. These are due 
to the relatively small box size of the 25 h' 1 Mpc run, which 
does not contain the massive groups and clusters in which 
gas can reach temperatures high enough to produce the high- 
est ionisation states of atoms with high atomic numbers. 
Variations for other ions are within 0.2 dex or less. 

A3 Resolution in the continuum spectrum 

In this Appendix we investigate how the resolution in wave- 
length used to interpolate the continuum tables affects the 
estimate of the contribution of the continuum to the total 
energy budget. 

CLOUDY provides the intensity of the continuum 
emissivity as a function of wavelength in the range A =0.01 
A-l m. The binning in wavelength is irregular, with smaller 
bins at shorter wavelengths. At full resolution, the contin- 
uum emissivity contains AW ~ 3700 elements. As shown in 
Fig. [4l the spectral shape of the continuum emissivity con- 
tains sharp edges at low temperatures. These need to be re- 
solved when interpolating the continuum to create full emis- 
sion spectra. However, for our purposes this might not be 
necessary, since we are only interested in estimating the total 
energy carried by the continuum. To speed up our calcula- 
tions, we have binned the continuum in a smaller number 



© 2013 RAS, MNRAS OOQ.mf22l 



22 Bertone et al. 




CP 

o 



40.2 - 



40.0 I i i i i i i i i i i i i i i i i i 

500 1000 1500 

number of bins in A 

Figure A3. Dependence of the specific emission (emission per 
unit volume) of the continuum on the number of bins used to 
integrate the intensity of the continuum. Differences in the results 
are well within 0.4 dex. 

of wavelengths than those contained in the tables, and we 
perform the interpolation over this new set of wavelengths. 

To verify that the binning does not affect the estimate 
of the specific emission of the continuum, we have tested 
how this varies as a function of AW- Results for runs with 
25 h' 1 Mpc at z = 2 and with 100 hT x Mpc at z = are 
shown in Fig.lA~3lfor N hin =50, 100, 200, 400, 800 and 1600. 
We find that a high number of bins does not substantially 
improve the accuracy in the estimate of the specific emission. 
Variations in the results are within 0.3 dex, and decrease 
to within 0.1 dex or less for Abm > 50. For this work, we 
assume A^in ~ 200 everywhere, except for Fig. [8] where we 
use AW = 100. 



APPENDIX B: LINE LISTS 

In this Appendix we provide the full lists of the individual 
emission lines that contribute most of the cooling energy. 
The fractions of the mean energy emitted per unit volume 
for individual lines as a function of temperature are shown 
in Figs. II ll and 1121 Only lines that contribute at least 1 per 
cent of the total emitted energy in at least one temperature 
bin are listed. Table [Bl] lists line for H, He, C, N, O, Ne and 
Mg. Table [B2] lists lines for Si, S, Ca and Fe. 

This paper has been typeset from a TpX/ I^TpX file prepared 
by the author. 
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Table Bl. List of the individual emission lines that account for most of the total emitted energy. The columns contain: i) the ionic 
species; ii) the wavelength of the line in A; iii) the temperature T m at which the line carries the largest fraction of the total energy; iv) 
and viii) the total fraction F of energy radiated through the line at z = and z = 2; v) and ix) the rank of the line according to the 
total fraction of radiated energy, vi) and x) the fraction F m of the energy carried by the line at the temperature T m where it carries the 
largest fraction of the total energy, at z = and z = 2; vii) and xi) the rank of the line according to the fraction of radiated energy at 
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Table B2. As Fig.[5T]- continued 
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